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This  is  the  final  report  on  a  project  for  investigating  the  dynamics  of 
aerosol  formation  in  gaseous  atmospheres  ranging  from  conditions  at  the 
Earth's  surface  to  those  in  the  extreme  upper  atmosphere.  The  research  in¬ 
volves  theoretical  and  experimental  determination  of  the  behavior  of  ultrafine 
aerosol  particles  at  high  concentrations.  The  normal  diffusion-limited  coagu¬ 
lation  and  growth  of  aerosols  is  expected  to  be  strongly  modified  as  the  par¬ 
ticle  diameter  and  the  interparticle  separation  approach  the  mean  free  path  of 
the  gas  molecules. 

In  the  first  year,  the  experimental  facility  was  set  up,  using  an  explod¬ 
ing  wire  generator  for  the  production  of  high  concentrations  of  metallic  aero¬ 
sols.  The  experiments  were  devoted  to  development  of  sampling  methods  for  the 
reduced  atmosphere  environment,  and  then  to  the  observation  of  the  types  of 
behavior  exhibited  by  single  component  aerosols.  The  types  of  aerosols  gener¬ 
ated  were  highly  irregular  in  appearance,  consisting  of  branched  chains  of 
small,  spherical  particles.  The  irregular  appearance  led  us  to  consider  the 
application  of  fractal  mathematics  for  describing  the  particles  and  the  physi¬ 
cal  processes  behind  their  growth. 

Much  of  the  second  year's  efforts  were  devoted  to  the  fractal  analysis  of 
the  aerosols  and  the  interpretation  of  the  results.  Fractal  analysis  is  still 
in  an  embryonic  state  with  regard  to  applications  in  physical  measurements, 
and  the  techniques  used  in  one  situation  may  not  be  the  best  to  use  in 
another.  The  investigation  of  techniques  for  measuring  fractal  dimensions  led 
us  to  reject  several  methods  as  inapplicable  to  the  problem  of  characterizing 
aerosol  particles  and  to  refine  the  use  of  others  to  obtain  consistency  among 
methods.  The  correlation  and  circular  dilation  techniques  have  been  the  meth¬ 
ods  of  choice  in  analyzing  micrographs  of  the  aerosol  particles.  These  meth¬ 
ods  also  belong  to  the  same  category  of  analytical  techniques  used  in  theoret¬ 
ical  models  of  particle  growth  processes,  which  allows  us  to  compare  the 
experimental  measurements  with  model  results. 
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In  the  third  year's  work,  emphasis  was  placed  on  answering  the  questions 
raised  in  the  experiments  and  analysis  of  the  prior  work.  Theoretical  test 
particles  were  developed  to  challenge  the  analytical  techniques.  Ensembles  of 
real  particles  were  measured  to  obtain  the  statistical  properties  of  the  frac¬ 
tal  dimensions  for  aerosols  grown  under  identical  conditions.  Several  new 
materials  were  used  for  the  aerosol  source  material  to  determine  what  differ¬ 
ences  of  growth  could  be  expected  and  what  other  effects,  such  as  electrical 
charging,  would  be  encountered. 

This  report  summarizes  the  work  of  the  prior  years  and  explains  in  detail 
the  work  carried  on  in  this  last  year's  effort. 


2.0  SUMMARY  OF  PRIOR  WORK 


2.1  YEAR  ONE 

Most  of  the  first  year's  effort  was  devoted  to  setting  up  the  experi¬ 
mental  apparatus  for  generating  and  growing  aerosols  at  reduced  pressures. 
While  this  was  being  done,  theoretical  studies  were  undertaken  of  the  kinds  of 
growth  that  could  be  expected. 

2.1.1  Theoreti cal  Work 

The  condensation  of  a  saturated  vapor  into  a  liquid  phase  of  droplets 
produces  a  quite  monodisperse  aerosol  when  the  growth  is  diffusion-limited. 
Subsequent  growth  of  the  aerosol  as  a  liquid  follows  general  patterns.  In  the 
continuum  regime,  particle  diameters  »  the  gas  mean  free  path,  the  particle 
size  distribution  evolves  toward  a  steady  shape  that  is  fairly  broad.  In  the 
free-molecular  regime,  particle  diameters  «  the  gas  mean  free  path,  the  size 
distribution  evolves  toward  a  similar  but  distinct  shape.  As  a  result  of  a 
study  of  the  growth  of  liquid  droplets  in  the  transition  regime,  between  the 
free-molecular  and  continuum  regimes,  a  third  size  distribution  was  determined 
that  is  significantly  narrower  than  the  other  two,  as  shown  in  Figure  1. 

This  transition  distribution  is  a  transient  distribution  that  will  even¬ 
tually  grow  into  the  continuum  regime  distribution;  however,  for  a  wide  range 
of  initial  particle  sizes  and  concentrations,  the  size  distribution  changes  so 
slowly  with  time  that  it  can  be  considered  quasistatic. 

When  the  liquid  aerosol  is  a  molten  metal  that  eventually  solidifies,  the 
frozen  droplets  may  then  continue  to  collide  and  stick  to  one  another,  but 
they  will  maintain  the  shape  and  diameter  they  had  at  the  time  of  freezing. 
This  fact  makes  an  analysis  of  the  diameters  feasible,  even  after  the  primary 
particles  have  been  incorporated  into  quite  complex  particle  agglomerates. 

This  will  be  taken  up  in  the  discussion  of  the  third  year's  work. 

Although  condensation  and  growth  of  aerosols  in  the  liquid  state  were 
given  a  theoretical  treatment  in  the  first  year  [1],  the  observations  of  the 
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experimental  particles  quickly  showed  that  the  solid  particle  growth  was  much 
more  complex  than  could  be  treated  with  the  standard  theoretical  methods. 

2.1.2  Experimental  Work 

The  aerosol  was  produced  with  an  exploding  wire  generator  (EWG)  [1]  on 
which  was  mounted  a  0.1  m^  growth  chamber  capable  of  being  evacuated  to  10  to 
4  atm  or  below.  The  aerosol  was  removed  from  the  chamber  for  capture  on  a 
filter  or  analysis  in  an  optical  particle  counter  (OPC)  by  capturing  it  in  a 
small  (145  cm^)  internal  chamber  that  was  then  raised  to  atmospheric  pressure 
and  flushed. 

The  chamber  technique  worked  satisfactorily  at  all  pressures.  At  low 
pressures,  the  aerosol  particles  settled  completely  in  about  0.5  h,  but  at 
1.0  atm,  large  numbers  of  particles  could  be  obtained  3  h  after  the  generating 
explosion.  A  shuttered  sampler  was  built  onto  the  floor  of  the  chamber  to 
capture  settling  particles  at  various  times  during  the  growth  period,  but  it 
was  rarely  used. 

The  OPC,  a  Particle  Measurement  Systems  LAS-X  model,  served  mainly  to 
monitor  the  particle  concentration.  The  counts  on  the  OPC  channels  served  as 
a  good  indicator  of  the  concentration  on  the  filters.  Pieces  of  the  filters 
were  mounted  on  stubs  for  observation  in  a  scanning  electron  microscope  iSEM), 
and  the  photographs  taken  on  the  SEM  became  the  primary  particle  measurement 
technique. 

The  experimental  results  in  the  first  year  were  a  series  of  particle 
pictures  taken  at  pressures  of  1.0,  0.1,  and  0.01  atm.  Attempts  were  made  to 
produce  particles  at  0.001  atm  without  success,  and  the  question  of  the  fail¬ 
ure  to  produce  particles  has  been  deferred  until  the  present  year.  The  par¬ 
ticle  pictures  that  were  taken  show  strikingly  irregular  structures,  the  kinds 
of  structures  that  defy  easy  description.  These  early  results  changed  the 
course  of  the  investigation  because  it  became  clear  that,  although  the  struc¬ 
tures  were  difficult  to  describe,  there  were  definite  pressure-related  differ¬ 
ences. 

Liquid  effects  are  important  during  the  first  few  milliseconds  of  growth, 
and  for  longer  times  at  lower  pressures,  but  such  effects  are  not  relevant  for 
most  of  the  aerosol  growth  we  have  observed.  Instead,  we  have  concentrated  on 
learning  how  to  explain  the  unusual  aspects  of  solid  aerosol  growth  and  apply¬ 
ing  those  lessons  to  the  changes  we  observed  with  pressure. 
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2.2  YEAR  TWO 

The  work  in  year  two  was  devoted  mainly  to  developing  an  understanding  of 
fractals  and  devising  computer  routines  to  assist  in  the  analysis  of  the  aero¬ 
sol  micrographs. 

2.2.1  Theoretical  Work 
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Fractals  are  self-similar  geometrical  figures  that  are  scale  invariant; 
they  have  the  same  appearance  no  matter  what  the  scale  of  viewing  is.  It  is 
clear  that  applying  fractal  analysis  to  real  aerosol  particles  involves  making 
certain  distinctions  in  the  fractal  terminology.  Exact  fractals  are  geometri¬ 
cal  structures  that  can  be  generated  by  an  algorithm,  and  natural  fractals  are 
structures  generated  by  some  random  process  that  follows  a  general  rule. 

Exact  fractals  can  conceptually  be  extended  to  all  scales  of  observation,  but 
natural  fractals  are  bounded  above  by  the  finite  size  of  the  structure  and 
below  by  the  fundamental  units  making  up  the  object.  Thus,  the  natural  frac¬ 
tal  structures  of  the  aerosols  have  a  lower  size  limit  imposed  by  the  primary 
spherical  particles  that  have  formed  the  aerosol  structure. 

The  natural  fractals  are  also  stochastic;  i.e.,  they  have  random  statis¬ 
tical  variations.  In  a  set  of  aerosols,  it  is  to  be  expected  that  individual 
members  will  have  somewhat  different  fractal  dimensions,  even  after  being 
generated  by  the  same  process.  The  average  of  the  individual  fractal  dimen¬ 
sions  can  be  expected  to  describe  the  characteristic  of  the  aerosol  struc¬ 
tures  . 

The  fractal  dimensions  associated  with  fractal  objects  is  a  measure  of 
the  structure  of  the  objects.  The  fractal  techniques  we  have  used  are  one¬ 
dimensional  descriptions  of  two-  and  three-dimensional  objects.  The  fractal 
dimensions  obtained  for  these  one-dimensional  descriptors  are  greater  than  1 
and  less  than  2. 

One  major  interest  in  using  fractals  to  characterize  these  aerosols  is 
that  considerable  modeling  work  has  been  done  on  the  fractal  dimensions  that 
result  from  different  growth  processes.  These  models  have  been  computer 
simulations  of  structures  growing  by  the  accretion  of  small  particles  onto  the 
larger  ones  or  by  the  aggregation  of  clusters  of  similar  size  with  one 
another,  under  the  influence  of  Brownian  (random  walk)  trajectories  or 
straight  line  trajectories. 
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The  models  clearly  separate  particle-cluster  and  cluster-cluster  growth 
according  to  fractal  dimension.  In  particular,  using  the  mass  fractal  defini¬ 
tion,  particles  that  collide  with  already  existing  clusters  by  Brownian  dif¬ 
fusion  exhibit  a  fractal  dimension  of  1.67  (5/3  theoretically).  In  higher 
dimensions,  this  is  extended  to: 

D  =  5/6d  ,  (2-1) 

where 

d  =  spatial  dimension  from  2  to  6  [3,4]. 

The  Brownian  collision  of  clusters  with  clusters  generally  leads  to  a 
fractal  dimension  of  1.45  (for  d  =  2)  or  1.75  (for  d  =  3)  under  a  wide  variety 
of  assumptions,  but  there  does  not  seem  to  be  an  exact  value.  These  charac¬ 
teristic  fractal  dimensions  also  are  generated  by  ballistic  (straight  line) 
trajectories,  mainly  because  clusters  cannot  interpenetrate  one  another  to  any 
great  extent. 

If  particle-cluster  collisions  occur  ballistically  (straight  line  trajec¬ 
tories),  the  fractal  dimension  is  about  1.95,  with  a  theoretical  limit  of  2 
for  infinitely  small  particles.  The  difference  between  1.95  and  2  occurs 
because  of  the  voids  between  tightly  packed  particles. 

When  dealing  with  images  of  particles  (strictly  two-dimensional),  the 
question  arises  of  how  to  relate  the  results  to  real  three-dimensional  par¬ 
ticles.  One  answer  is  that  for  sparse  particles,  in  which  there  is  relatively 
little  screening  of  parts  of  the  image  by  other  parts,  the  fractal  dimension 
should  be  the  same  in  two  or  in  three  dimensions  [5].  Another  answer  is  that 
the  fractal  dimension  of  the  particle  must  obey  a  relation  derived  from  cau¬ 
sality  [6]: 

Do  >  d  +  1  -  Di  ,  (2-2) 

where 

Dq  =  particle  fractal  dimension 

Dj  =  fractal  dimension  of  the  trajectory  of  an  incoming  particle. 

If  the  trajectory  is  Brownian,  then  Dj  is  equal  to  the  space  dimension.  If  it 
is  ballistic,  the  D]  is  1.  Therefore,  the  particle  fractal  dimension  must 
fall  between  the  rather  broad  limits  of  1  and  2  for  two-dimensional  processes 
or  between  1  and  3  for  three-dimensional  processes. 
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2.1.2  Experimental  Work 


In  our  use  of  fractals,  we  started  with  a  tiling  method,  described  in 
Reference  [1],  but  had  difficulty  interpreting  the  results  of  an  analysis  of 
exact  fractals.  A  tile  could  be  counted  as  occupied  either  if  a  large  frac¬ 
tion  of  it  contained  image  pixels,  or  a  very  small  fraction.  The  resulting 
fractal  dimensions  showed  large  variations,  depending  on  the  occupancy  rule. 
Tiling  is  a  process  of  approximating  the  area  of  an  image,  and  other  methods 
approximate  the  linear  dimensions  of  the  image.  Tiling  may  prove  useful  in 
future  investigations  of  aerosol  particles  in  which  area  coverage  is  impor¬ 
tant,  but  it  was  abandoned  in  this  one. 

The  remaining  development  of  the  analysis  methods  concentrated  on  three 
broad  categories:  perimeter  measurements,  dilation  measurements,  and  correla¬ 
tion  measurements.  These  have  been  recognized  techniques  for  computing  frac¬ 
tals,  but  there  has  been  little  standardization  of  the  methods  for  analyzing 
images . 

2.2.2. 1  Perimeter  Method-- 

In  this  method,  the  perimeter  of  an  object  is  estimated  by  counting  the 
number  of  steps  of  a  given  size  required  to  circumscribe  the  object,  and  plot¬ 
ting  it  against  the  step  size  on  a  log-log  plot.  The  fractal  dimension  of  the 
perimeter  is  l+lm|,  where  m  is  the  slope  of  the  plot,  assuming  that  there  is  a 
sufficiently  long  linear  region  to  say  that  the  object  is  indeed  fractal  in 
nature.  An  equation  for  expressing  the  perimeter  as  a  function  of  the  step 
size  is: 

P(r)  ~  rl-D  ,  (2-3) 

where 

D  =  fractal  dimension,  between  1  and  2  in  two  dimensions. 

The  perimeter  (or  external  hull)  fractal  does  not  take  into  account  an 
object  with  interior  perimeters,  such  as  we  see  frequently  in  our  micrographs. 
Furthermore,  the  perimeter  fractal  is  not  easily  suited  to  computer  analysis. 
On  the  other  hand,  the  perimeter  fractal  does  seem  to  be  useful  in  picking  out 
changes  of  shape  between  the  primary  spheres  making  up  an  object  and  the  gen¬ 
eral  shape  of  the  object  itself.  At  small  scales,  the  perimeter  approaches 
the  geometric  perimeter  of  the  circles  making  up  the  particle  image,  and  the 


fractal  dimension  is  close  to  1 .  At  larger  scales,  the  perimeter  bridges  gaps 
around  the  particle,  and  the  structure  changes  more  rapidly  with  scale,  giving 
a  larger  fractal  dimension.  We  have  observed  some  particles  that  were  much 
longer  along  one  axis  than  the  other,  and  these  particles  show  a  return  to  a 
fractal  dimension  of  1  at  scales  between  the  length  of  the  minor  axis  and  the 
length  of  the  major  axis.  Thus,  in  that  range  the  particle  behaves  as  if  it 
were  a  one-dimensional  structure. 

2. 2. 2. 2  Dilation  Method-- 

A  commonly  used  alternative  definition  of  the  fractal  dimension  is  based 
on  the  distribution  of  mass  around  the  center  of  mass  of  the  object.  If  the 
object  is  fractal,  the  mass  is  given  by: 

M(r)  ~  rD  ,  (2-4) 

where 

D  =  fractal  dimension  (less  than  the  spatial  dimension). 

We  have  called  this  a  mass  fractal,  but  more  accepted  usage  is  the  term  "den¬ 
sity  fractal."  The  method  used  to  calculate  it  is  the  dilation  method  because 
a  series  of  boxes  is  expanded  around  the  center  of  mass  to  determine  the  mass 
within  each  box. 

In  this  case,  r  represents  a  "radius"  of  gyration  around  the  center  of 
mass.  This  fractal  dimension  definition  is  quite  intuitive  in  its  application 
because  it  corresponds  closely  with  the  fraction  of  area  obscured  by  a 

particle.  A  sparsely  distributed  particle  tends  to  have  a  low  fractal  dimen¬ 

sion,  closer  to  1  than  2,  and  a  densely  distributed  particle  has  a  high  frac¬ 
tal  dimension,  approaching  2. 

It  has  been  accepted  practice  to  use  square  boxes  in  the  method,  but  we 
think  that,  for  comparison  with  other  analysis  methods,  the  expansions  should 
be  with  circles.  The  difference  between  using  boxes  and  circles  is  that  the 
box  method  confuses  the  scale  r  by  a  factor  of  1.4,  blending  some  parts  of  the 

image  that  are  farther  from  the  center  than  others. 

The  expansion  around  the  center  of  mass  means  that  all  the  small-scale 
information  is  obtained  from  the  region  near  the  center  of  mass.  This  ne¬ 
glects  small-scale  information  in  other  parts  of  the  particle,  unless  auxil¬ 
iary  expansions  are  performed  over  all  the  particle. 


At  large  scales,  approaching  the  boundaries  of  the  particle,  the  dilation 
method  approaches  a  constant  mass,  the  total  mass  of  the  particle,  and  the 
fractal  dimension  no  longer  has  meaning.  This  boundary-limited  behavior  oc¬ 
curs  at  smaller  scales  when  expansions  are  performed  around  centers  other  than 
the  center  of  mass.  Boundary  effects  on  the  mass  fractal  must  be  anticipated 
in  the  analysis,  and  the  boundary-limited  regions  must  not  be  given  a  fractal 
interpretation. 

2. 2. 2. 3  Correlation  Method-- 

The  correlation  function  is  a  relation  describing  the  relative  frequency 
of  finding  two  parts  of  the  image  at  a  given  distance  from  one  another.  It  is 
an  alternative  way  of  calculating  the  density  fractal.  For  fractal  objects, 
the  correlation  function  is  related  to  the  fractal  dimension  by: 

C(r)  ~  r^-d  ^  (2-5) 

where 

D  =  fractal  dimension 

d  =  spatial  dimension. 

Because  the  fractal  dimension  is  less  than  the  spatial  dimension,  the  correla¬ 
tion  function  decreases  with  r. 

Programs  were  written  to  calculate  the  correlation  function  directly  and 
involved  scanning  the  image,  multiplying  pixel  values,  and  summing  at  various 
values  for  r.  The  process  is  manageable  for  small  images  and  small  values  of 
r,  but  the  time  consumed  in  the  calculation  grows  quickly  for  larger  images. 
The  correlation  method  does  give  a  very  complete  picture  of  the  fractal  char¬ 
acteristics  of  the  particles  at  small  scales  and  seems  to  agree  well  with  the 
dilation  method  at  intermediate  scales. 

Both  the  mass  and  correlation  fractals  show  similar  behavior  for  a  given 
particle  and  range  of  scales.  The  correlation  function  itself  decreases  with 
r,  but  if  it  is  multiplied  by  r^,  then  the  result  is  very  similar  to  the  mass 
fractal : 

r2.c{r)  ~  rD-drd  =  rd  ,  (2-6) 

because  d  =  2  for  the  images.  In  this  way,  it  is  feasible  to  use  the  dilation 
method  to  obtain  the  fractal  characteristics  at  large  scales,  where  it  is  a 
fairly  efficient  technique,  and  the  correlation  method  at  small  scales,  where 
it  is  effective. 


3.0  SUMMARY  OF  PUBLISHED  WORK  FROM  THE  CURRENT  YEAR 


The  fractal  analysis  methods  as  described  have  been  adequate  for  the 
initial  investigation.  The  dilation  analysis  has  been  expanded  to  include  the 
circle  dilation  concept,  where  circles  are  the  expansion  element,  rather  than 
boxes.  Several  particles  of  various  shapes  have  been  compared  with  analyses 
by  circle  dilation,  box  dilation,  and  correlation.  The  results  are  being 
prepared  for  publication,  as  noted  on  the  list  of  publications  [6]  in 
Section  7.0. 

In  general,  the  overall  agreement  between  the  box  and  circle  dilations  is 
good.  However,  as  the  dilations  approach  the  size  of  the  particles,  the  two 
methods  give  different  results.  Generally,  when  there  is  a  change  in  the 
slope  of  the  log-log  curve,  the  change  shows  up  more  distinctly  and  at  a  smal¬ 
ler  scale  in  the  circle  dilation  method  than  in  the  box  method.  We  attribute 
this  to  the  more  clearly  defined  scale  lengths  of  the  circle  dilation. 

The  circle  dilation  and  the  correlation  methods  also  generally  agree  over 
most  of  their  common  scale  range,  and  the  changes  of  slope  occur  at  approxi¬ 
mately  the  same  scale  lengths.  For  some  particles,  the  correlation  and  dila¬ 
tion  methods  give  indistinguishable  results,  but  for  others  grown  under  the 
same  conditions,  the  agreement  is  obviously  poorer. 

The  sensitivity  of  the  fractal  analysis  methods  to  changes  in  the  shapes 
of  the  particles  was  the  subject  of  substantial  investigation.  This  consisted 
of  a  three-part  program:  analysis  of  multiple  particles  grown  under  identical 
conditions;  analysis  of  simulated  three-dimensional  particles  under  different 
projections;  and  analysis  of  simulated  three-dimensional  mass  distributions. 

The  density  (mass)  fractals  of  two  groups  of  particles  were  measured  on 
single  micrographs  taken  at  pressures  of  1.0  and  0.1  atm.  The  micrographs 
were  at  a  relatively  low  magnification,  but  the  images  were  enlarged  optically 
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with  only  a  slight  degradation  of  the  resolution  that  could  be  attained  with 
high  magnification  micrographs.  The  methods  used  were  the  correlation  and 
circle  dilation.  The  following  table  (Table  1)  summarizes  the  results. 

All  of  the  resulting  fractal  dimensions  are  equal  within  the  standard 
deviations  of  the  averages.  These  numbers  are  consistent  with  the  cluster- 
cluster  growth  in  two-dimensions  mentioned  earlier.  The  visual  appearance  of 
these  particles  is  also  suggestive  of  clusters  or  chains  of  particles  coming 
together.  There  are  a  few  particles  in  each  group  whose  fractal  dimensions 
are  quite  low  or  quite  high  for  cluster-cluster  aggregation.  The  appearance 
of  these  deviant  particles  is  not  particularly  unusual  in  comparison  with  the 
rest.  The  results  of  this  analysis  confirm  the  use  of  the  density  fractal  as 
a  consistent  measure  of  the  particle  shape. 

The  second  test  involved  the  creation  of  an  artificial  particle.  This 
could  have  been  done  with  a  physical  model,  but  instead  was  done  with  a  com¬ 
puter  algorithm.  The  algorithm  computed  the  centers  of  primary  particles  in 
three-dimensional  space  in  a  way  that  mostly  retained  the  center  of  mass  near 
the  axes'  origin,  but  it  gave  a  reasonably  random  character  to  the  particle. 
The  particle  could  be  viewed  under  different  orientations  and  analyses  in 
those  projections.  Under  many  orientations,  the  artificial  particle  resembled 
natural  particles  whose  fractal  dimensions  had  been  previously  determined.  In 
those  cases,  the  density  fractal  had  a  very  similar  dimension  to  the  natural 
particles. 

The  size  of  the  primary  spheres  making  up  the  artificial  particle  could 
be  varied  at  will.  Values  were  chosen  to  (1)  ensure  that  adjacent  positions 
in  the  particle  overlapped  substantially;  (2)  cause  occasional  gaps  between 
the  primary  spheres;  or  (3)  to  cause  the  primary  spheres  to  be  mostly  isolated 
from  one  another.  These  cases  are  shown  in  Figure  2,  where  the  primary  sphere 
size  varies  across  the  columns  and  the  orientation  varies  down  the  columns. 

The  numbers  at  the  bottoms  of  the  columns  show  the  mean  and  standard  deviation 
of  the  dimensions  in  each  column  determined  by  the  circle  dilation  and  corre¬ 
lation  methods. 

Again,  the  two  methods  are  in  generally  good  agreement  with  one  another. 
However,  there  is  a  strong  variation  of  the  fractal  dimensions  as  the  size  of 
the  primary  particles  decreases.  This  variation  is  indicative  of  the  nature 
of  the  fractal  analysis,  namely  that  the  object  should  be  scale  invariant  over 


TABLE  1.  FRACTAL  DIMENSIONS  FOR  PARTICLE  ENSEMBLES 


Pressure 

(atm)  Number  of  particles  Dd 
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Figure  2.  An  artificial  particle  shown  with  different  orientations  and 
primary  particle  sizes. 


a  range  of  scales.  The  artificial  introduction  of  gaps  between  the  particles 
breaks  the  scale  invariance  at  a  scale  comparable  to  the  gap  length.  This  can 
even  be  detected  in  some  of  the  plots  of  the  correlation  and  dilation  func¬ 
tions  . 

The  orientation  effects  with  this  artificial  particle  are  within  the 
range  of  variation  of  the  natural  particles  and  may  be  assumed  to  be  negligi¬ 
ble.  The  question  of  orientation  will  be  addressed  with  natural  particles  in 
Section  4.0. 

The  last  test  of  the  sensitivity  of  the  fractal  analysis  involved  another 
set  of  algorithm-generated  mass  distributions,  but  this  time  with  no  attempt 
to  construct  contiguous  spherical  masses.  Instead,  mass  distributions  were 
generated  in  three  dimensions  by  distributing  point  masses  in  a  uniformly 
random  fashion  over  the  surfaces  of  concentric  shells.  The  number  of  point 
masses  in  each  shell  was  determined  by  a  power-law  relation  between  the  total 
mass  and  the  shell  radius.  The  power-law  exponents  varied  between  0.1  and  20 
for  the  distributions.  The  point  of  the  test  was  to  see  the  ability  of  the 
fractal  analysis  to  recover  the  three-dimensional  mass  distribution  from  the 
two-dimensional  projected  image. 

The  test  cases  and  correlation/di latfon  plots  are  shown  in  Figure  3  for 
eight  different  power-law  relations.  This  test  should  be  counted  a  failure 
because  neither  the  correlation  nor  the  dilation  could  reproduce  the  generat¬ 
ing  power-law  exponent.  The  correlation  function  was  almost  invariant  for  all 
the  power  laws  used,  and  no  information  could  be  obtained  from  it.  The  dila¬ 
tion  function,  on  the  other  hand,  shows  behavior  that  can  be  used  to  distin¬ 
guish  between  different  power  laws.  The  fractal  dimension  computed  from  the 
slope  of  the  curve  does  not,  however,  relate  even  semiquanti tati vely  to  the 
generating  exponents. 

These  results  show  the  main  difference  between  the  dilation  and  correla¬ 
tion  methods.  The  correlation  method  examines  all  points  in  the  space  around 
a  given  point,  while  the  dilation  method  examines  all  points  at  a  fixed  dis¬ 
tance  from  the  dilation  center.  The  correlation  method  will  pick  up  the  in¬ 
fluence  of  the  gaps  between  mass  points  much  more  effectively  than  can  the 
dilation  method,  whereas  the  dilation  method  will  discern  the  overall  mass 
distribution  with  radius.  When  the  mass  points  occupy  contiguous  space,  the 
two  methods  give  very  similar  results. 


In  light  of  Mafidelbrot '  s  statement  about  the  fractal  dimensional  i  t>  of  a 
projection  of  a  fractal  [6],  it  is  not  surprising  to  find  that  power-la^  ex¬ 
ponents  greater  than  2  give  strange  results  in  the  projection.  The  dilation 
plots  for  generating  exponents  less  than  2  do  give  extended  ranges  of  pcwer- 
law  behavior,  consistent  vvith  fractal  character.  However,  the  fractal  dimen¬ 
sion  does  not  reflect  the  generating  power-law  exponent. 

Further  reflections  on  the  character  of  three-dimensional  fractals  led  us 
to  conclude  that  all  projections  of  such  objects  will  have  dimensions  of  about 
2  when  the  three-dimensional  structure  has  a  dimension  greater  than  2. 


4.0  DETAILS  OF  UNPUBIISHED  WORK  OF  THE  CURRENT  YEAR 


The  remainder  of  the  final  year's  effort  has  been  devoted  to  answering 
some  of  the  questions  that  have  occurred  in  the  course  of  the  work.  Broadly 
speaking,  these  are  the  questions  of  the  universality  of  the  phenomena  we  have 
observed,  the  problems  of  generating  particles  at  low  pressures,  the  effects 
of  electric  fields  on  the  growth  of  the  aerosol,  the  effects  of  covapor i zat i on 
of  other  materials  with  the  metal,  the  nature  of  the  fractal  analysis  under 
actual  rotation  of  the  specimens,  and  the  primary  particle  size  distribution 
under  various  gas  pressures. 

We  will  consider  these  questions  and  the  answers  produced  in  the  experi¬ 
ments  in  the  following  sections. 

4.1  CONFIRMATORY  RESULTS  WITH  OTHER  METALS 

Our  approach  for  confirming  the  universality  of  the  results  we  obtained 
with  silver  aerosol  was  to  use  a  variety  of  other  metals  under  the  same  condi¬ 
tions.  This  meant  using  metals  with  higher  and  lower  melting  boiling  points, 
which  would  simultaneously  allow  us  to  explore  electrical  charging  of  the 
aerosols  during  the  explosion  phase  and  the  low  pressure  formation  problems. 

Time  limitations  prevented  a  total  survey  of  the  metals  that  were  avail¬ 
able.  Table  2  shows  the  ones  tested  and  the  associated  physical  properties 
that  are  important. 

This  selection  vy.ers  a  Aide  range  of  melting  and  bciling  points.  A 
previous  report  [7]  had  pointed  out  our  difficulties  in  generating  aerosols 
with  aluminum.  These  attempts  had  been  made  in  air  atnr:  spheres ,  and  'is  pa'- 
ticles  were  observed,  either  visual  1>  or  with  an  QpC.  .-.hen  the  a  uminum 
wastried  once  more  in  a  nitrogen  atmosphere,  the  aerosol  formed  v,as  comparable 
visually  to  that  formed  from  other  metals,  strongly  suggesting  that  the  pre¬ 
vious  failures  were  due  to  the  oxidation  of  the  molten  or  vap'tiz!=>d  a  i  umi ’iliT:. 


TABLE  2.  METALS  USED  FOR 

AEROSOL  GENERATION 

Element 

Melting  point 
(“C) 

Boiling  point 
(“C) 

Densi ty 
(g.'  cm3) 

Si  1 ver 

961 

1,950 

10.5 

Indi urn 

156 

1,450 

7.3 

T  i  n 

232 

2,260 

7.3 

Zinc 

419 

907 

7.1 

Molybdenum 

2,620 

4,507 

10.2 

P1 ati num 

1,773 

4,300 

21.5 

I  ron* 

1,535 

3,000 

7.9 

Al umi num* 

660 

2,057 

2.7 

*Aerosols  from 

these  metals  were  not  col 

llected,  but  simply 

observed 

visual ly . 


With  the  exception  of  the  zinc,  the  metals  all  formed  spherical  pi  Miiary 
particles.  The  zinc  fucmed  crystallites  SGmc/.'hat  larger  in  mean  size  than  tlie 
others.  The  subsequent  coagulation  of  all  the  materials  was  similar.  Mit'':- 
graphs  of  the  various  metals  under  different  conditions  are  shown  in  Appen¬ 
dix  B . 

The  density  fractal  dimensions  of  these  materials  are  summarized  in 
Table  3. 

4.2  ANALYSIS  OF  LOW-PRESSURE  PROBLEMS 

After  repeatedly  failing  to  generate  and  observe  any  particles  at  all 
with  silver  at  0.001  atm,  we  raised  the  pressure  until  particles  were  observed 
at  0.002  atm.  Other  metals  exhibited  similar  behavior,  with  complete  Failure 
to  generate  particles  below  0.002  atm.  There  are  individual  differences  from 
metal  to  metal,  but  we  have  localized  the  critical  pressure  to  within  only  a 
factor  of  2 . 

Understanding  the  explosion  process  itself  may  help  identify  the  cause 
for  this  behavior.  The  details  of  the  explosion  follow  the  general  steps 
1  i sted  below  [3] : 

1.  The  flow  of  current  heats  the  wire  to  its  melting  point.  Inertia 
and  magnetic  fields  contain  the  molten  wire  and  allow  current  to 
continue  flowing. 

2.  The  wire  is  heated  well  above  its  vaporization  point.  It  is  super¬ 
heated  because  the  time  scale  is  too  short  for  equilibrium  boiling. 
In  a  reduced  gas  pressure,  it  may  not  heat  to  the  same  temperature 
as  in  a  high-pressure  atmosphere. 

3.  Vaoor  pockets  form  in  the  liquid,  but  heating  continues  as  Icna  as 
tnere  is  a  continuous  path  in  the  liquio. 

4.  Vaporization  disrupts  the  current  flow.  Hie  remaining  liquid  I’ia- 
ments  coalesce  into  spherical  droplets. 

5.  Thermionic  emission  causes  the  dr-plets  to  emit  electrons  and  become 
positively  charged.  Because  the  vapor  density  is  high,  electee' 
multiplication  is  suppressed. 

6.  The  vapor  expands  rapidly.  The  liquid  droplets  are  carried  intc' 
regions  of  decreasing  vapor  concentration.  In  the  decreasing  vapor, 
electron  multiplication  takes  place,  and  the  arc  may  re-strike. 


TABLE  3.  DENSITY  FRACTALS  FOR  METALLIC  AEROSOLS 


Element 

Pressure 

(atm) 

Dc 

Dd 

Si  1 ver 

1.0 

1 

.40  ±  0.15 

1.36  ±  0.15 

0.1 

1 

.35  t  0.06 

1.37  ±  0.09 

0.01 

1 

.70  i  0.02 

1.60  ±  0.02 

0.002  — 

Minimum 

for  formation 

Indium 

1.0 

1 

.38  ±  0.08 

1.41  ±  0.13 

0.01 

1 

.47  ±  0.20 

1.42  ±  0.24 

0.005 

1 

.54  ±  0.12 

1.54  ±  0.11 

0.002 

1 

.49  ±  0.11 

1.47  ±  0.10 

0.005  -- 

Minimum 

for  significant 

formation 

Tin 

1.0 

1 

.58  ±  0.12 

1.47  ±  0.12 

0.01 

1 

.51  ±  0.06 

1.49  ±  0.08 

0.002  -- 

Minimum 

for  formation, 

only  single 

particles 

Zinc 

1.0 

1.38 

1.45 

0.01 

1.59 

1.51 

0.004  — 

Minimum 

for  formation 

Molybdenum 

1.0 

1.56 

1.53 

0.01* 

1.68 

1.74 

0.002 

1.29 

1.33 

0.002  - 

Minimum 

for  formation 

Platinum 

1.0 

1.40 

1.55 

0.01* 

1.56 

1.64 

0.002 

1.30 

1.40 

0.002  -- 

Minimum 

for  formation 

'See  the  section  on  low-pressure  problems  (Section  4.2)  for  additional 
informat i on . 


This  analysis  of  the  wire  breakup  has  been  supported  by  high-speed  X-ray 
photographs  of  wires.  The  mechanisms  leading  to  breakup  may  not  be  unique, 
but  it  is  clear  that  some  liquid  droplets  remain  from  the  original  wire.  It 
is  not  totally  vaporized  by  the  electrical  current. 

Now,  the  superheated  liquid  droplets  expand  into  the  sample  chamber  where 
competing  processes  take  place.  First,  there  is  cooling  by  the  ambient  gas. 
There  is  potential  continued  evaporation  of  the  material  from  the  droplets, 
and  there  is  potential  condensation  of  metallic  vapor  onto  the  droplets  or 
spontaneous  nucleation  and  growth  of  new  droplets. 

The  minimum  radius  that  can  support  continuous  growth  is  given  by  the 
Kelvin  equation: 

rc  =  2Ms, RTdlnS  ,  (4-1) 

where 

r^  -  particle  radius 

M  =  molecular  weight  of  metal 

s  =  surface  tension  of  metal  (850  dyne 'em) 

R  =  gas  constant 
T  =  temperature  of  metal 
d  =  density  of  metal 
S  =  saturation  ratio  [9], 

If  the  radius  is  less  than  this  critical  value,  the  droplet  will  ulti¬ 
mately  evaporate.  If  the  radius  is  larger,  it  will  grow,  'he  proolem  in 
applying  this  relation  is  the  great  uncertaintr  about  the  experimental  situa¬ 
tion. 

Under  the  steps  leading  to  the  formation  given  above,  the  temperature  c' 
the  liquid  drop  must  be  at  least  equal  to  the  normal  boiling  point  of  sii.e-", 
2,483  K.  If  the  superheating  is  substantial,  the  temperature  might  be  as  high 
as  2,600  K.  It  is  also  not  clear  that  the  ambient  pressure  could  affect  this 
temperature  strongly.  The  saturation  ratio  is  definitely  an  equilibrium  con¬ 
cept:  the  ratio  of  ambient  partial  pressure  to  the  vapor  pressure  at  the  sur¬ 
face  of  the  particle. 
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Calculations  of  the  critical  radius  over  a  range  of  temperatures  corre¬ 
sponding  to  boiling  at  0.001  atm  up  to  7  atm  il,560  to  5,040  K)  produce  a  very 
narrow  range  of  critical  sizes.  For  S  =  1.01,  the  critical  radius  decreases 
from  0.18  to  0.05  /im.  For  S  =  1.02,  the  range  is  0.09  to  0.03  /im,  and  for  S  = 
1.05,  the  range  is  0.04  to  0.01  /^m.  These  results  are  right  in  the  range  of 
primary  oarticles  encountered.  Any  particles  formed  with  smaller  radii  than 
these  will  evaporate.  Any  larger  ones  will  grow.  Table  4  summarizes  the 
critical  radius  (diameter)  calculations  for  a  number  of  materials  and  includes 
mean  primary  particle  diameters  measured  on  micrographs  of  the  particles. 

The  calculated  critical  diameters  and  measured  primary  particle  sizes  are 
in  reasonably  good  agreement  for  saturation  ratios  of  1.05  and  are  generally 
in  poor  agreement  for  saturation  ratios  of  1.01.  The  critical  diameters  for 
the  various  metals  are  reasonably  close  to  one  another  at  a  given  saturation 
ratio.  This  agreement  is  strengthened  in  the  experimental  measurements  of 
particle  diameters,  where  relatively  minor  differences  exist.  The  outstanding 
exception  is  for  zinc  at  1.0  atm.  The  measured  diameter  is  considerably 
higher  than  for  the  other  metals,  or  for  zinc  at  lower  pressures.  In  the 
micrographs,  the  primary  particles  of  zinc  appear  to  be  small  crystals  rather 
than  solidified  spheres.  This  crystal  1 ization  suggests  that  the  formation 
process  is  significantly  different  and  probably  involves  the  slow  cooling  of 
droplets  that  had  not  completely  vaporized. 

The  uncertain  variables  in  this  calculation  are  the  surface  tension, 
which  is  rarely  available  over  the  temperature  range  of  interest,  and  the 
equilibrium  vapor  pressure.  The  use  of  a  constant  saturation  ratio  implicitly 
assumes  that  the  ambient  pressure  sets  the  vapor  pressure  of  the  vaporized 
metal  within  the  ratio  factor.  Thus,  at  0.01  atm,  we  have  assumed  the  vapor 
pressure  of  the  metal  is  7.6  Terr  and  that  the  temperature  of  the  vaoor  is 
equal  to  the  equilibrium  temperature  at  that  pressure.  The  measurements  of 
primary  particle  diameter  are  also  susceptible  to  sampling  and  statistical 
error  because  the  micrographs  are  at  their  limit  of  resolution  and,  at  most, 
contain  a  few  hundred  primary  particles. 

Nonetheless,  the  calculations  show  that  the  primary  particles  tend  to  be 
larger  as  the  pressure  is  lowered,  and  the  saturation  ratio  of  1.05  is  of  the 
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TABLE  4.  CRITICAL  DIAMETERS  FOR  METALLIC  AEROSOL  FORMATION 

Element 

Pressure 

(atm) 

Surface 
tension 
(dyne  'em) 

Di ameter 
S=1.01 
ilim) 

Diameter 

S=1.05 

(/im) 

Di  ameter 
measured 
(/im) 

Si  1  ver 

1.0 

850 

0.17 

0.03 

0.060 

" 

0.1 

850 

0.20 

0.04 

0.045 

II 

0.01 

850 

0.24 

0.05 

0.061 

Indi urn 

1.0 

515 

0.23 

0.05 

0.053 

" 

0.01 

515 

0.25 

0.05 

0.064 

M 

0.004 

515 

0.26 

0.05 

0.050 

" 

0.002 

515 

0.26 

0.05 

0.067 

T  i  n 

1.0 

480 

0.17 

0.03 

0.069 

" 

0.01 

480 

0.22 

0.04 

0.042 

II 

0.002 

480 

0.24 

0.05 

0.038 

Zinc 

1.0 

763 

0.29 

0.06 

0. 126 

II 

0.01 

763 

0.40 

0.08 

0.049 

M 

0.004 

763 

0.043 

0.09 

0.041 

Molybdenum 

1.0 

2,080 

0.20 

0.04 

0.064 

11 

0.01 

2,080 

0.25 

0.05 

0.054 

II 

0.002 

2,080 

0.28 

0.06 

0.047 

Platinum 

1.0 

1,700 

0.16 

0.03 

0.040 

II 

0.01 

1,700 

0.22 

0.05 

0.049 

II 

0.002 

1,700 

0.25 

0.05 

0.047 

Al umi num 

1.0 

860 

0.04 

0.01 

Iron 

1.0 

1,800 

0.21 

0.04 

Copper 

1.0 

1,200 

0.14 

0.03 

Magnesi urn 

l.O 

502 

0.25 

0.05 

*Count  median  diameter  (CMD). 

right  magnitude.  From  these  results,  theoretical  estimates  of  primary  par¬ 
ticle  size  can  be  made,  as  has  been  done  in  Table  4  for  aluminum,  iron, 
copper,  and  magnesium. 

This  prediction  of  the  primary  particle  sizes  still  does  not  answer  the 
low-pressure  formation  problem.  Examination  of  the  chamber  walls  after  each 
explosion  at  0.001  atm  showed  the  presence  of  spots  of  the  metal  that  could 
only  have  been  deposited  in  the  molten  state.  Different  metals  showed  dif¬ 
ferent  diameter  spots.  At  the  slightly  higher  pressures  for  which  aerosol 
formation  occurred,  there  was  still  evidence  of  the  spots,  but  it  was  not  as 
pi  enti  fill  . 

The  conclusion  that  we  reach  from  these  consistent  observations  is  that 
the  exploded  metal  droplets  collide  with  the  walls  of  the  chamber  before  evap- 
porating.  The  flight  path  is  approximately  55  cm  from  the  wire  to  the  wall 
where  the  spots  were  observed.  Raising  the  gas  pressure  to  0.002  atm  stops 
most  of  the  droplets  for  most  metals  within  this  distance,  allowing  them  to 
evaporate.  Although  the  stopping  distance  does  depend  on  gas  and  particle 
density,  particle  radius,  and  initial  velocity,  it  is  worthwhile  to  examine 
the  relations  that  can  be  derived. 

The  stop  distance  for  high-velocity  particles  is  well  above  the  Stokes' 
law  range  and  is  given  by; 

s  =  4  3  fip'flgd  dRe  C[)Re  ,  (4-2) 

where 

s  =  distance 
Pp  =  particle  density 
fig  =  gas  density 
d  =  particle  diameter 

Re  =  particle  Reynolds  numbe'-,  which  depends  on  velocity 

Cq  =  drag  coefficient,  which  also  depends  on  Reynolds  number  [9]. 

The  Reynolds  number  for  particles  is  given  by: 

Re  =  vd^g  /i  ,  \4-3) 

where 

V  -  particle  velocity 
d  '  pa rt i c I e  d i ameter 
fi  ■  gas  V i scos i t \ . 


The  drag  coefficient,  Cq,  is  a  dimensionless  quantity  that  depends  on  the 
Reynolds  number.  For  the  intermediate  turbulent  range,  2  <  Re  <  800,  the  drag 
coefficient  varies  as  Re'^'^,  and  for  the  turbulent  range,  Re  >  1,000,  it  is 
constant  for  all  practical  purposes. 

Assuming  that  the  droplets  are  all  generated  v\ith  the  same  dimaeters  and 
initial  velocities,  the  stopping  distance  will  mainly  depend  on  the  inverse  of 
the  gas  density.  Therefore,  the  stop  distance  will  vary  inversely  with  pres¬ 
sure  . 

From  this  result,  a  change  of  pressure  by  a  factor  of  2  can  make  the 
difference  between  droplets  stopping  entirely  within  the  chamber  or  being 
deposited  on  the  chamber  walls.  If  we  assume  a  stop  distance  of  50  cm  at 
0.002  atm,  then  we  can  extrapolate  to  higher  and  lower  pressures.  The  results 
are  shown  in  Table  5. 

We  see  from  the  table  that  laboratory  production  of  aerosol  by  an  explod¬ 
ing  wire  generator  will  become  impractical  below  0.001  atm  because  the  size  of 
the  chamber  required  to  stop  the  droplets  in  the  gas  is  excessive.  The  rate 
of  coagulation  and  growth  of  the  primary  particles  is  also  a  function  of  their 
concentration,  once  formed.  As  the  gas  pressure  is  reduced,  even  assuming 
that  the  stopping  distance  is  within  the  confines  of  the  chamber,  the  average 
concentration  decreases  with  the  third  power  of  the  stopping  distance. 

One  final  problem  that  had  been  anticipated  originally  appeared  at  low 
pressures.  For  the  molybdenum  and  platinum  wires,  the  explosion  failed  to 
destroy  the  wire.  In  fact,  only  a  very  small  fraction  cf  the  mass  was  euipc- 
rated,  and  the  wires  showed  no  signs  of  meltinq.  Aerosol  pai  tides  were  ptc- 
duced,  although  at  much  lower  concentrations  than  expected. 

The  explanation  for  the  persistence  of  the  wires  is  that  the  surce  c' 
current  through  the  wires  caused  rapid  heatina,  but  the  resdtance  or  the 
wires  also  increased  with  temperature.  This  allowed  the  voltage  across  the 
wire  to  increase  to  the  point  that  gas  breakdown  occurred,  shunting  the  I'e- 
maining  energy  from  the  capacitor  into  the  gas.  The  breakdown  strength  of  the 
gas  is  proportional  to  gas  pressure,  and  this  breakdown  would  be  expected  to 
occur  first  at  the  lowest  pressures.  The  high  melting  points  of  the  platinum 
and  molybdenum  allow  the  temperatures  and  resistances  to  increase  substantial¬ 
ly  without  losing  tne  inteqritv  of  the  wire. 
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TABLE  5.  bTOP  DISTANCE  AS  A  TUNtFION 
OF  PRESSURE 


Pressure  Stop  distdiice 

(atm)  (cm) 


10.0 

0.01 

1.0 

0.1 

0.1 

1 

0.01 

10 

0.001 

100 

0.0001 

10-5 

10-5 

1,000 

10^ 

105 

27 


4.3  MAJOR  ELECTRIC  FIELD  LFFECTS 


Metallic  aerosols  are  particularly  susceptible  to  electric  field  effects 
because  even  electrically  neutral  particles  can  experience  strong  dipole 
polarization  and  torques  in  an  electric  field.  The  concentration  of  an  elec¬ 
tric  field  at  the  extremities  of  a  conducting  particle  causes  strong  gradient 
forces  there,  which  leads  us  to  expect  that  the  random  agglomeration  patterns 
will  become  much  more  linear.  This  has  been  borne  out  in  computer  models  of 
agglomeration  with  long-range  forces  [10]. 

Initial  observations  of  different  materials  in  the  EWG  using  a  Plexiglas 
chamber  showed  that  some  materials,  particularly  iron,  formed  extensive 
chains,  while  others,  such  as  silver  or  copper,  did  not  [1].  The  suspected 
low  charge  on  the  silver  aerosol  was  confirmed  with  Faraday  cage  measurements 

The  visual  observation  of  chain  formation  is  considered  a  good  indicator 
of  electric  fields  influencing  aerosol  growth.  Flowever,  in  the  metal  chamber 
the  observed  chain  growth  was  confined  to  the  metal  posts  supporting  the  wire 
before  the  explosion  and  the  bolt  heads  nearby.  The  probable  reason  for  this 
phenomenon  is  as  follows. 

Any  cloud  of  charged  aerosol  particles  is  self-dispersive.  That  is,  the 
charges  on  the  particles  causes  them  to  repel  one  another  and  the  cloud  ex¬ 
pands.  The  equation  governing  the  motion  is  quite  general  and  easy  to  calcu¬ 
late  if  the  particle  mobilities  are  known.  In  the  Plexiglas  chamber,  the 
charged  particles  are  quickly  deposited  on  the  chamber  walls,  but  they  do  not 
lose  their  charge.  From  those  positions  on  the  walls,  their  charges  create 
e'ectric  fields  in  the  chamber  that  continue  to  influence  the  growth  of  the 
I'emaining  neutral  particles  and  agglomerates .  In  the  metal  chamber,  on  the 
other  hand,  the  particles  caught  on  the  walls  lose  their  charges  instantlv  an 
can  no  longer  influence  ttie  aggl  omerat  i  rn . 

The  electrodes  and  the  bolts  near  them  can  have  residual  potential^  on 
them,  due  to  the  circuitrv,  and  so  aie  in  a  position  to  collect  part'.’es  .n 
dendritic  chains,  if  enough  charge  exists  initially.  The  presence  of  den¬ 
dritic  growths,  and  their  degree,  is  also  considered  a  measure  of  the  charge 
on  the  aerosol  particles  from  the  explosion.  Table  6  summarizes  the  observa¬ 
tions  of  dendrite  growth. 

The  observation  of  charge  effects  is  rr.ughlv  in  proportion  to  the  meltin 
points  and  not  to  tiie  boiling  pi.ints.  Even  v.’th  f'tp  variation  c  *  irot  t  o'' i  - 1  -  , 


TABLE  6.  METALLIC  DENDRITE  OBSERVATIONS 


Melting  point 
(“C) 

Boi 1 i ng  poi nt 
(°C) 

Dendr i te 
formation 

Si  1 V  e  r 

961 

1,950 

Very  si i ght 

Indi urn 

156 

1,450 

None 

Tin 

232 

2,260 

None 

Z  i  nc 

419 

907 

Very  slight 

Molybdenum 

2,620 

4,507 

Moderate 

Platinum 

1  ,  773 

4,300 

Moderate 

I  ron 

1,535 

3,000 

Considerable 

A1 umi num 

660 

2,057 

None 

Note:  Approximately  five  times  the  normal  amount  of  material  xas  e-.ploded 
for  the  iron  wire,  leading  to  the  heavier  dendrite  formation. 


none  of  the  cent i meter- 1 ong  chains,  which  .were  so  obvious  in  the  Plexiglas 
chamber,  was  observed  in  the  metal  chamber.  The  longest  particle  that  was 
observed  was  roughly  1  mm  in  length  (for  molybdenum).  Therefore,  an  addi¬ 
tional  electric  field  source  was  placed  in  the  chamber  through  a  Plexiglas 
plate  for  observation. 

The  source  was  a  1,8  in.  diameter  rod  connected  to  a  high  voltage  power 
supply.  When  this  rod  is  centered  in  the  12-in.  diameter  chamber,  the 
electric  field  at  any  radial  distance  from  the  center  is  given  by; 

£(r)  =  V/[r  ln(rc  rr)]  - 

where 

r(-  =  chamber  radius 

r^  =  rod  radius. 

Voltages  of  200  to  20,000  V  were  applied  to  the  rod  after  a  molybdenum  wire 
explosion.  The  application  of  any  voltage  within  this  range  led  to  dramatic 
formation  of  chains  of  up  to  a  centimeter  in  length,  but  additional  effects 
were  observed. 

The  chains  either  carried  an  initial  charge  or  formed  in  contact  with 
either  the  chamber  wall  or  the  rod  because  there  was  a  steady,  rapid  drift  of 
chains  toward  or  away  from  the  rod.  At  the  higher  applied  voltages,  the 
chains  would  accelerate  in  one  direction  or  another,  but  they  would  be 
arrested  in  flight  as  a  piece  of  the  chain  broke  off  and  carried  the  attrac¬ 
ting  charge  to  the  electrode.  The  remainder  of  the  chain,  carrving  a  repul¬ 
sing  charge,  was  driven  in  the  opposite  direction.  This  motion  of  the  chains 
had  the  effect  of  scavenging  the  chamber  of  small  particles  ve-'v  quicklv,  ana 
the  concentrat ion  dropped  many  times  faster  than  it  would  have  done  withr,.: 
the  action  of  the  electric  field. 

In  th^s  case,  the  electric  field  was  strong  enough  for  the  differentia' 
forces  on  the  ends  of  the  chains  to  cause  the  breakage  of  them.  Tiie  constant 
electric  field  also  caused  an>  particles  in  contact  with  the  wall  or  t'od  tc* 
acquire  an  induced  charge.  For  many  particles,  the  force  on  the  induced 
charge  was  strong  enough  to  lift  them  from  the  surface  into  flight,  but  the 
majority  remained  on  the  walls  pe>'manent  1  y . 

These  observations  strcuiglv  support  the  role  of  electric  fields  in  tl'e 
formation  of  the  chains.  Ihe  sitsjnqer  the  field,  the  taster'  tne  ,-hain- 


both  because  of  the  attraction  of  neutral  particles  to  the  end  of  the  chain 
and  because  of  the  motion  of  the  chain  through  the  relatively  stationary  r.eu- 
tral  cloud. 

Although  time  was  not  available  for  us  to  collect  aerosol  particles  grown 
under  these  conditions,  the  agglomerated  aerosols  from  the  higher-melting- 
point  materials  do  show  evidence  of  more  chainlike  growth  than  the  low-melt¬ 
ing-point  materials  do. 

4.4  COVAPORIZATION  CONFINEMENT 

Covaporization  of  a  metal  and  another  material  was  a  part  of  the  pro¬ 
jected  work  in  this  last  year's  effort.  The  reason  for  doing  this  is  that 
many  aerosol -generating  schemes  always  involve  more  than  one  component,  either 
deliberately  or  through  contamination.  The  shock  from  the  exploding  wire 
would  have  enough  energy  to  vaporize  or  shatter  a  foil  wrapped  around  the 
wire,  and  it  would  produce  a  mixed  aerosol  environment  in  which  the  gi'cwth  of 
the  metallic  aerosol  could  be  substantially  modified. 

An  experimental  explosion  was  performed  with  a  silver  wire  inserted  into 
a  small  piece  of  Teflon  tubing  before  being  mounted  in  the  generator.  Al¬ 
though  the  Teflon  was  chosen  mainly  for  convenience,  a  dielectric  such  as 
Teflon  would  have  a  likelihood  of  acquiring  charges  during  disruption  and 
affecting  the  subsequent  growth  of  the  aerosol. 

The  explosion  of  the  wire  did  vaporize  much  of  the  tubing,  but  several 
pieces  were  recovered  from  the  chamber  floor  after  the  explosion.  Low-power 
microscopic  examination  of  the  larger  pieces  showed  that  molten  droplets  i-.e’"e 
embedded  in  them,  much  like  the  molten  droplets  that  had  impacted  tne  :!'ambe'' 
walls.  Small  particles  of  material,  presumably  Teflon,  ..e'-e  obse’.ed  t.:  he 
attracted  to  the  electrodes  after  the  explosion,  tut  the  m.inbe’'  did  '  t 
iticrease  after  the  first  few  minutes. 

In  all  other  respects,  the  growth  of  the  siKer  aeroscl  h;,’ij,.ed  cte 
usual  course,  with  the  particles  growing  in  size  until  they  could  te  dist’:’- 
guished  visually,  and  then  decreasing  slowly  by  settling.  There  was  no  macrc- 
scopic  evidence  of  altered  growth.  Time  prevented  the  pursuit  of  microscopic 
samples. 

4.5  FRACTAL  STABILITY  UNDER  ROTATION 

The  problem  of  using  fractal  analysis  methods  -jv.  three-d  i  mens  i  :i'a  1  y-i- 
tides  in  two-d i mens i onal  micrographs  was  aopr-ached  in  Iwi  wav'.  The  ti'  : 
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is  reported  in  Reference  [7],  where  an  artificial  particle  was  computer-gener¬ 
ated  in  a  way  that  allowed  it  to  be  rotated  and  projected  from  any  point  of 
view.  Several  views  of  this  particle  are  shown  in  Figure  3. 

Although  the  apparent  fractal  dimension  changes  with  the  size  of  the 
primary  particles,  viewing  across  the  rows,  the  fractal  dimensions  under  rota¬ 
tion,  viewing  down  the  columns,  are  in  reasonable  agreement.  The  standard 
deviation  of  the  different  orientations  is  comparable  to  the  standard  devia¬ 
tion  for  particle  ensembles  in  Table  1. 

The  second  approach  was  to  view  some  naturally  formed  particles  at  dif¬ 
ferent  orientations  (indium  at  1.0  atm).  This  was  done  by  tilting  the  SEM 
stage  at  various  angles  before  taking  micrographs.  (This  technique  is  often 
used  to  make  stereoscopic  pairs  for  the  detection  of  three-dimensional  fea¬ 
tures.)  A  group  of  indium  particles  was  found  that  offered  a  variety  of  par¬ 
ticle  shapes;  pictures  were  taken  with  the  electron  beam  at  angles  of  0,  30, 
50,  and  74  degrees  with  respect  to  the  normal.  The  fractal  dimensions  for 
five  of  these  particles  are  presented  in  Table  7.  Two  of  the  particles  are 
shown  in  Figure  4,  where  it  can  be  seen  that  the  character  of  the  particles 
does  not  change  strongly  for  most  orientations.  The  scales  of  these  figures 
are  not  all  the  same,  as  they  would  be  if  taken  directly  from  the  micrographs. 
These  are  edited  images  taken  from  the  computer  screen. 

The  results  of  this  test  are  a  good  indicator  of  the  error  in  the  density 
fractal  dimension  that  could  be  expected  from  orientation  effects.  It  is  of 
the  same  order  as  the  expected  variation  from  particle  to  particle.  The 
standard  deviations  compare  favorably  with  those  for  the  theoretical  particle 
orientation. 

We  conclude  that  the  analysis  of  two-dimensional  images  can  be  used  suc¬ 
cessfully  to  obtain  the  three-dimensional  fractal  dimensions  of  aerosol  par¬ 
ticles,  at  least  as  long  as  the  fractal  dimension  is  strictly  less  than  2. 

4.6  PRIMARY  PARTICLE  SIZE  DISTRIBUTIONS 

Visual  inspection  of  the  micrographs  gives  the  impression  that  the  pri¬ 
mary  particle  size  becomes  smaller  and  more  varied  as  the  pressure  is  lowered. 
During  the  digitization  of  the  primary  particles,  the  size  distributions  were 
computed.  A  simple  measure  of  the  width  of  the  distribution  is  given  by  the 
ratio  of  the  mass  median  diameter  (MMD)  to  the  count  median  diameter  (CMD) . 

The  MMD  is  computed  as  the  diameter  that  divides  the  total  particulate  mass 


TABLE  1. 

ORIENTATION 

EFFECTS 

FOR  NATURAL 

PARTICLES 

Part  id  e 

n 

0 

30 

50 

74 

Avg 

1 

Fc: 

1.39 

1 .44 

1.43 

1.41 

1.42  ±  0.03 

Fq: 

1.33 

1 .40 

1.40 

1.41 

1.38  ±  0.05 

2 

Fc: 

1.30 

1.26 

1.32 

1.31 

1.30  ±  0.03 

Fd-- 

1.33 

1.55 

1.49 

-- 

1.46  ±  0.14 

3 

Fc: 

1.46 

1 .49 

1.48 

1.58 

1.50  ±  0.06 

Fd: 

1.49 

1.55 

1.55 

1.49 

1.52  ±  0.04 

4 

Fc: 

1.31 

1.18 

1.29 

1.45 

1.31  ±  0.13 

Fd: 

1.40 

1.20 

1.41 

1.71 

1.43  ±  0.24 

5 

Fc: 

1.79 

1.76 

1.79 

1.78  ±  0.02 

Fd: 

-- 

1.85 

1 .81 

1.96 

1.87  ±  0.10 

R.M.S.  deviation  0.11 


R.M.S.  =  Root  mean  square. 
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into  two  equal  halves,  and  the  CMD  is  the  diameter  that  divides  the  particles 
into  two  groups  of  equal  numbers. 

If  the  primary  particles  are  all  the  same  size,  then  the  HMD  and  the  CMD 
are  equal.  As  the  distribution  includes  more  particle  sizes,  the  MMD  grcws 
faster  than  does  the  CMD  because  of  the  volume  dependence  of  mass.  Therefore 
the  more  polydisperse  the  distribution  is,  the  greater  the  ratio  between  MMD 
and  CMD.  Table  8  summarizes  the  results  for  all  the  metals  measured. 

From  the  table,  it  cm  be  seen  that  the  silver  and  indium  results  are 
responsible  for  the  impressions  cited  above.  Both  have  a  much  more  nearly 
monodisperse  primary  distribution  at  1.0  atm  than  the  others,  and  the  ratio 
increases  at  the  lower  pressures.  The  other  metals  give  a  much  less  clear 
i ndi cation . 

It  is  tempting  to  speculate  that  the  initial  size  distribution  is  nearly 
monodisperse  at  all  pressures  and  that  the  increasing  polydi  specs  i  ty  at  low-er 
pressures  is  the  result  of  liquid  (molten)  particle  coagulations  in  the  first 
few  seconds  of  growth  before  solidification.  The  coagulation  process  should 
lead  to  more  polydisperse  primary  particles,  ultimately  growing  into  the  dis¬ 
tributions  shown  in  Figure  1.  Although  we  think  this  interpretation  is 
qualitatively  correct,  there  is  not  enough  evidence  to  conclude  so 
quanti tati vely . 

4.7  DISCUSSION 

Before  coming  to  conclusions  about  the  results  presented,  we  would  like 
to  summarize  the  computer  model  results  publisned  for  two-  and  three-dimen¬ 
sional  aggregation  (Tables  9  ana  10,  respective.,,)- 

Now  that  we  are  confident  that  the  two-dimens''onal  fractal  analysis  c'.e 
us  the  fractal  dimensions  in  three  dimensions,  fc r  D  <  C,  we  must  abandc' 
interpretations  of  the  aerosol  behavior  based  on  two-dimensional  models.  The 
two-dimensional  models  describe  the  aggregation  of  particles  confined  tv;  o 
plane,  such  as  on  the  surface  of  a  liquid.  Their  use  is  i nappropr i ate  for 
aerosol  formation  with  three-dimensional  growth  characteristics.  This  means 
that  the  range  of  fractals  that  we  have  encountered,  with  dimensions  from 
about  1.2  to  1.8,  must  be  measured  against  the  results  for  the  three-dimen¬ 
sional  models. 

Concentrating  on  the  Table  !0,  we  see  that  only  the  c 1  us ter-c 1  us  ter 
aggregation  results  give  fractals  in  the  proper  range  to  describe  the  fra.vta' 


lABLE  9.  ^RAC^A1.  DIMENSIONS  OBIAINl-D  t  ROM 
TWO-DIMENSIONAL  AGGREGAIION  MODE  IS 

Model  D 


Linear  trajectory, 

particle-cluster  1.95  ±  0.002 

Brownian  trajectory, 

particle-cluster  (theory)  1.6667 

Brownian  trajectory, 

particle-cluster  1.73  ±  0.06 

Linear  trajectory, 

cluster-cluster  1.60  ±  0.05 


Brownian  trajectory 
cluster-cluster 


1.44  ±  0.02 


TABLE  10.  LRACTAL  DIMENSIONS  OBTAINED  FROM 
THREE-DIMENSIONAL  AGGREGATION  MODELS 


I 


Model 


Linear  trajectory, 
particle-cluster 


Brownian  trajectory, 

particle-cluster  (theory) 


Brownian  trajectory, 
pa  rt i c 1 e-c 1  us  ter 


Linear  trajectory, 
cl uster-cl uster 


2.97  ±  0.08 


2.5000 


2.45  ±  0.1 


1.81  t  0.03 


P 


m 


v.V'C'l 
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we  have  eiicouiitered.  The  nature  of  cluster-cluster  aggregation  obscures  the 
differences  between  the  trajectories  because  clusters  tend  to  penetrate  each 
other  very  little,  whether  or  not  they  travel  in  straight  lines. 

Our  computer  analysis  results  for  distributions  of  particles  show  that 
very  few  differences  can  be  distinguished  for  power-law  mass  distributions 
with  exponents  of  2  or  greater.  This  is  consistent  with  the  lemma  that  frac¬ 
tals  of  dimensionality  <2  can  be  measured  in  two  dimensions.  During  the 
course  of  the  experiments,  we  have  encountered  clumps  of  particles  that  prob¬ 
ably  result  from  gravitational  scavenging  of  small  particles  by  large  clus¬ 
ters.  These  have  a  dense  appearance  and  a  density  fractal  dimension  near  2 
(1.^3  by  correlation,  1.84  by  dilation.)  Presuming  that  these  clumps  exhibit 
a  three-dimensional  fractal  dimension  of  about  2.97  (for  linear  trajectory, 
particle-cluster  interactions),  the  fractal  dimension  of  their  projection  is 
indeed  almost  equal  to  2. 

This  means  that  the  fractal  analysis  of  images  cannot  distinguish  any 
three-dimensional  characteristics  with  dimensionality  greater  than  2,  i.e. 
particle-cluster  aggregation.  The  conclusion  we  must  then  draw  is  that  frac¬ 
tal  dimensions  of  1.6  to  1.8  in  the  micrographs  represent  cluster-cluster  ag¬ 
gregation  only.  This  is  reasonable  in  light  of  the  predominance  of  clusters 
in  the  micrographs,  but  it  leaves  unanswered  the  mechanism  responsible  for  the 
measured  fractals  in  the  range  1.2  to  1.5  that  are  generally  more  prevalent 
than  are  the  higher  values. 

The  lower  fractal  dimensions  indicate  that  some  mechanisms  are  present 
that  are  not  accounted  for  in  the  standard  models.  The  apparent  fractal  d’- 
mension  in  the  range  of  two-dimensional  cluster-cluster  aggregation  causes  us 
to  consider  the  possibility  of  forces  that  confine  the  motions  of  the  par¬ 
ticles  to  two  dimensions.  These  might  be  either  electrical  forces  or  sheai 
forces . 

Shear  forces  arise  when  a  velocity  gradient  exists  in  the  gas.  The  part 
of  the  particle  in  the  higher  velocity  region  attempts  to  move  faster  ti'an  the 
part  in  the  lower  velocity  region.  This  results  in  an  orientation  force  that 
tends  to  move  the  particles  into  the  higher  velocity  region.  Consider  two 
particles  moving  in  the  gas  near  one  another.  If  the  velocity  shear  is  essen¬ 
tially  planar,  as  occurs  for  gas  moving  near  a  flat  wall,  and  the  particles 
are  at  different  distances  f  I'om  the  wall,  they  will  travel  with  s  i  gn  i  t  i  ca;' 1 1  . 
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different  velocities  and  spend  little  time  near  each  other.  Therefore,  the 
cha  :es  of  their  diffusing  to  one  another  and  sticking  are  small.  If, 
however,  they  are  at  the  same  distance  from  the  wall,  they  will  spetid  a  long 
time  near  each  other,  even  though  both  are  moving  with  respect  to  the  wall, 
and  have  a  much  higher  probability  of  colliding. 

If  the  particles  are  moving  at  slightly  different  distances  from  the  wall 
and  slightly  different  velocities,  then  they  may  collide  in  a  grazing  con¬ 
figuration  and  stick  together.  This  also  has  the  effect  of  confining  the 
relative  motion  to  a  plane  because  any  cross-flow  movements  would  reduce  the 
probability  of  collision  greatly. 

Thus,  motion  in  a  shear  flow  can  have  the  effect  of  confining  the  dif¬ 
fusive  motion  to  a  plane  and  might  result  in  a  large  fraction  of  tne  particles 
having  the  two-dimensional  cluster-cluster  characteristic. 

Electric  fields  have  been  shown  to  have  a  pronounced  effect  on  the  par¬ 
ticles  in  causing  the  chain  structures  to  grow.  Electric  fields  also  play  a 
role  in  orienting  complex  particles.  The  presence  of  an  electric  field 
induces  charge  separation  (polarization)  in  any  metallic  conductor.  Positive 
and  negative  charges  move  apart  until  their  separation  causes  a  field  that 
exactly  cancels  the  external  field  in  the  interior  of  the  conductor. 

These  separated  charges  can  be  acted  on  by  the  external  field  to  create 
torques  on  the  particle.  The  torques  tend  to  orient  the  particle  so  that  the 
longest  axis  of  the  particle  points  in  the  direction  of  the  field.  Any  dif¬ 
fusive  motion  between  particles  will  tend  to  take  place  along  the  field  cirec- 
tion  in  preference  to  the  other  directions  because  the  pol  ari  zat  i  on  charce-;-  w' 
one  particle  tend  to  attract  the  polarization  charges  on  anothe-'  particle 
along  the  field  lines.  If  the  particles  approach  each  other  on  parai’ei  T-tr'.: 
lines,  the  polarization  chai'qes  tend  to  repel  each  other,  "fiiv  is  the  rFO'.  r 
for  the  formation  of  chains,  but  it  also  tends  to  bring  the  t.vo  oarticles  'rt' 
contact  at  their  extremities. 

Of  the  two  forces,  we  think  the  electric  field  mechanism  is  more  likely 
because  the  evidence  of  particle  charges  is  already  confirmed  by  the  formation 
of  dendrites  and  chains.  Any  shear  motion  would  have  to  be  driven  by  the 
turbulence  immediately  after  the  explosion  or  by  convective  motions  near  the 
walls  after  the  gas  in  the  chamber  has  become  quiescent. 
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still,  the  appcrtt  rince  of  the  low  fractal  diiiiension  in  the  particles  pro¬ 
vides  an  indication  of  the  yiowth  processes  that  occlh  red  during  the  lifetime 
of  the  particle. 

4.3  CCMCLUSIONS 

The  conclusions  we  can  draw  from  this  work  are  as  follows: 

1.  Metallic  aerosol  can  be  formed  from  a  molten  or  vapor  state  in  the 
laboratory  down  to  a  pressure  of  about  0.001  atm.  Formation  at 
lov.er  pressures  becomes  exceedingly  difficult  and  may  require  much 
larger  active  volumes  and  initial  mass  of  material. 

2.  Metallic  aerosol  particles  can  be  produced  in  an  inert  gas  atmos¬ 
phere  from  the  molten  state.  In  the  presence  of  o>ygen,  the  aer'sol 
will  generally  be  an  oxide  (nonmetal  1 ic) . 

3.  The  melting  point  of  the  material  controls  the  degree  of  charges 
associated  with  the  disruptive  process  in  the  e'>ploding  wire  gener¬ 
ator.  In  another  form  of  generator,  such  as  a  thermal  source,  the 
charging  might  not  be  as  significant,  but  it  would  be  expected  to 
follow  the  same  order  for  materials. 

4.  The  primary  particles  formed  are  usually  spherical  and  are  consis¬ 
tent  with  a  condensational  growth.  The  average  diameter  of  the 
particles  is  in  good  agreement  with  theory,  using  a  saturation  ratio 
of  1.05,  approximately.  There  is  surprisingly  little  variation  of 
primary  particle  size  from  metal  to  metal,  but  exceptions  may  be 
found,  namely,  aluminum. 

5.  There  is  a  slight  trend  tcv.ard  forming  smaller  primary  particles  as 
the  gas  pressure  is  lowered.  There  is  a  definite  trend  toward  form¬ 
ing  more  polydisperse  primary  particles  as  the  gas  pressure  is 
lowered.  This  may  arise  from  standard  processes  involving  liquia 
droplet  coagulation  in  which  tiie  prima-'v  particles  -  ema  i  n  nr:Uen  f  >■ 
longer  per-iods  at  reduced  pressures. 


Fractal  analysis  methods  for  micr'ogt'aphs  of  aerosol  particles  a<e 
now  wel  ' -developed.  Correlation  and  dilation  techi'icues  cjive  re¬ 
sults  in  good  agreement  with  one  anotner.  The  perimeter  t'-arta' 
atialvsis  is  less  appropriate  for  aerosol  gt'owth  stucie^. 

Fractal  analvsis  of  particles  does  aistinguish  differences  in  crc-.t'' 
mechanisms.  Without  the  guidance  of  models,  however,  the  infor'ma- 
tion  obtained  would  be  difficult  to  relate  to  growth  processes. 

The  fractal  analysis  of  images  gives  reliable  fractal  dimensions  for 
characteristic  fractals  less  than  2  in  three  dimensions.  It  cannot 
distinguish  fractal  characteristics  that  would  have  dimensions 
greater  tlnni  ?. 
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Fractal  analysis  is  not  sensitive  to  the  orientation  of  the  par¬ 
ticles,  within  reason.  The  variation  of  fractal  dimension  i.ith 
orientation  is  of  the  same  oider  as  from  particle  to  particle  within 
an  ensemble  of  particles. 

For  the  aerosols  studied,  three-dimensional  cluster-cluster  aggrega¬ 
tion  (D  -  1.7)  is  the  dominant  growth  mechanism. 

There  is  substantial  evidence  of  microscopic  charge  effects  in  the 
aerosol  growth,  represented  by  chain  formation  and  low  fractal  di¬ 
mensions  (D  -  1.4). 

Gas  pressure  affects  the  density  fractal  dimension.  It  generally 
increases  as  the  pressure  is  reduced.  This  means  that  the  aerusol 
structures  are  more  compact  at  lower  pressures.  The  trend  is  toward 
ballistic  aggregation  at  pressures  below  0.001  atm,  but  this  has  not 
been  verified  experimentally. 


5.0  IMPLICATIONS  OF  THIS  WORK  FOR  SPACECRAFT  SURVIVABILITY 


In  this  section,  we  discuss  some  of  the  implications  of  this  program  on 
the  overall  spacecraft  survivability  issue. 

The  generation  and  growth  of  metallic  aerosols  in  a  vacuum  environment 
would  require  one  of  two  conditions:  (1)  a  generation  technique  substantially 
different  from  the  exploding-wi re  generator,  or  (2)  an  artificial,  low-pres¬ 
sure  atmosphere.  The  experiments  have  shown  that  a  critical  saturation  ratio 
is  needed  to  foster  the  growth  of  primary  particles  and  that  a  finite  time  is 
required  for  the  molten  droplets  produced  in  the  EWG  to  evaporate  and  cool. 

This  time  is  available  in  the  laboratory  because  the  droplets  are  halted 
by  the  ambient  gas.  In  the  orbital  vacuum,  the  cooling  would  be  by  radiation 
only,  and  it  would  take  a  significantly  longer  time.  During  that  cooling 
time,  the  particle  cloud  would  have  expanded  far  beyond  the  dimensions  we  ha.e 
used  in  the  laboratory. 

A  different  generator  principle  would  allow  the  dispersion  of  primary 
particles  into  the  vacuum,  but  their  subsequent  growth  would  still  be  oeter- 
mined  alotig  the  lines  of  this  investigation.  Aggregation  of  primary  pat  tides 
in  a  vacuum  would  be  by  ballistic  trajectories,  without  question.  Theref^M'e, 
the  aggregates  would  be  compact,  almost  solid  bodies,  unless  some  extertv.d 
forces  cause  some  preferential  collision  directions. 

Orbital  mechanics  v.'ould  '.reate  a  shearing  velocity  field:  Those 
ticles  in  lower  orbits  v.  i  11  move  with  a  higher  velocity,  and  these  in  r'gne'' 
orbits  will  move  with  a  lov.er  velocity.  Grazing  contacts  between  part'des 
would  give  a  tendency  toward  chains  and  clusters,  but  with  a  much  lower  rate 
than  diffusive  motion  would. 

The  elongated,  open  structures  of  the  aerosol  particles  we  have  produced 
would  be  effective  electromagnetic  radiation  scatterering  centers  because 
their  cross-section  is  roughly  equivalent  to  a  solid  of  the  same  diameter  at 
the  wavelengths  of  interest.  These  structures  cannot  form  without  some  d'^'- 
fusive  motion  unless  the  clusters  making  them  up  are  pt'eformed .  L\en  it  das- 


ter  subunits  were  dispersed  into  the  vacuum,  growth  would  be  slow  without 
diffusion  or  differential  velocities  to  cause  collisions. 

After  exposure  in  the  vacuum  environment  for  a  period  of  time,  the  par¬ 
ticles  would  begin  to  be  charged  by  the  ionized  species  present  and  by  photo¬ 
electric  emission  from  exposure  to  sunlight.  The  development  of  a  net  uni¬ 
polar  charge  would  hasten  the  dispersal  of  the  particulate  cloud,  whereas  the 
development  of  bipolar  charges  from  the  photoelectric  effect  would  encourage 
chain  formation.  Other  factors  influencing  the  rate  of  dispersal  of  the  cloud 
are  its  initial  ejection  velocity  from  the  parent  craft  and  the  difference  in 
orbital  velocities  between  the  parent  and  the  cloud,  and  between  parts  of  the 
cloud. 

The  maximum  control  of  the  growth  and  dispersal  of  a  particulate  cloud 
would  occur  in  a  confined  environment  with  a  suitable  gas  atmosphere.  The 
concept  of  a  large  balloon  is  the  most  obvious  method  of  maintaining  this 
control.  With  a  captive  gas  atmosphere,  the  boundaries  of  the  initial  dis¬ 
persal  can  be  set  by  the  stopping  distance  of  the  initial  particles.  The 
subsequent  coagulation  and  growth  rates  can  be  controlled  by  determining  the 
concentration  of  particles  and  gas  pressures.  Containment  of  the  cloud  also 
minimizes  problems  of  spurious  contamination  and  of  removing  the  cloud  when 
desi red . 
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